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Introduction 

The estimated errors of peak areas are primarily 
important results for the complete analytical pro- 
cess (e.g., activation analysis or environmental ra- 
diocontamination assay). Many excellent program 
routines were offered in the literature for determin- 
ing peak areas; however, no unambiguous solution 
has been given so far for assigning accurate error 
estimation methods. 

As an introduction to our "uniform concept" 
some basic principles are to be stated: 

- The error estimation would establish an accu- 
rate confidence limit around the calculated re- 
sults: this region would include the "true 
value." Thus the error is a measure of accuracy 
rather than that of precision. 

- The above condition is met only if the esti- 
mated error combines all statistical (random) 
and systematic (bias) errors. 

- The methods for peak area determination and 
error estimation, respectively, would be closely 
interdependent. 

- The actual formulation is strongly related to the 
given detection system. 
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is the total count rate of the multiplet de- 
termined by summing and baseline sub- 
traction, 

is the number of overlapping peaks 
is the count rate of the i'th peak, 
is the net peak height, and 
is the peak width. 



Peak heights are computed from least squares fit 
with fixed and previously calibrated widths. The 
error of the i'th count rate is given by eq (2): 
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The first two terms represent the "bias" error, the 
third one contains the "random" error. Peak height 
errors Aa are generated in the least squares fit. 

The count rate / for singlets is the maximum of 
the summed and baseline-subtracted net counts (/s) 
and the integral of the fitted peak shape function 
(/p). If we consider this step as a "one component 
apportioning," the error of the count rate is given 
by eq (3): 
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The resemblance of eqs (2) and (3) is quite clear. 

If the radioactivities attributed to the identified 
isotopes are calculated not individually but in a 
common least squares procedure (called interfer- 
ence correction), the estimated peak area errors 
must be applied as weight factors in the least 
squares solution. 



Theory 

Of course it is necessary that a uniform method 
be used for calculating peak areas of singlets and 
multiplets as well. The "apportioning" method 
proved satisfactory for resolving overlaps: 
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Results 

Numerous reference materials were analyzed to 
confirm the reliability of both the peak area 
computation and error estimation methods (such as 
Bowen's kale, NBS and IAEA reference materi- 
als). In testing the experimental standard devia- 
tions, it was revealed that some analytical results 
were biased. However, most results were still reli- 
able because the known true values in most cases 
did fall into the calculated confidence intervals. 
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Summary 



We elaborated a uniform formulation for calcu- 
lating estimated count rate errors in gamma spec- 
trometry. Experimental results confirmed that our 
method gave accurate confidence limits combining 
both random and bias errors. 
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1. Introduction 

Trace amounts of H, C, N and O can be deter- 
mined by only a few methods, often with rather 
poor accuracies. Charged particle activation analy- 
sis (CPAA) is highly reliable for C, N and O, and is 
used for calibration of other methods. Elastic recoil 
detection analysis (ERDA) and nuclear reaction 
analysis (NRA) have recently been utilized for H 
and D near the surface. Accuracies in these analy- 
ses are discussed on the basis of our experimental 
data. 



2. Accuracy in CPAA 

2.1 Advantage and Error in CPAA 

The advantage of CPAA for C, N and O consists 
of inherent high sensitivity, freedom from various 
contaminations, and reliability in the use of com- 
parators. The following reactions are used for acti- 
vation: ■=C('He,a)"C or '^C(d,n)"N; '*N(p,a)"C; 
and '*0(^He,p)'*'F. Uncertainty in the measurement 
of bombarding particle beam current frequently be- 
comes a major cause of error in CPAA. Also, some 
bombardment damage of the sample often results in 
noticeable overestimation of C and O. 



2.2 Precision for O and C 

We analyzed O in a silicon wafer repeatedly nine 
times from 1982 to 1985 by the nondestructive 



measurement of the '^F annihilation radiation with 
a Ge(Li) detector [1]. The results gave the mean 
value of 10.04 ppm wt with or =0.1 8 ppm. 

We examined chemical separation of "C in Si, 
using "C-containing Si prepared by proton bom- 
bardment of B-doped Si [1]. The following separa- 
tion method was selected: alkali-dissolution of the 
pulverized sample; KMnO* oxidation of the "C un- 
der microwave heating to 700 °C; and generation 
of "CO2 and its conversion into Li2"C03. The mea- 
surement of "C activity before and after the sepa- 
ration and of carrier recovery indicated that 
9.0+2.0% of the "C apparently disappeared. Our 
results obtained by this separation has thus been 
corrected for this 9%. 

For the separation of '^N, dry fusion into '■'N2 
and wet distillation as ''NHs are used [1]. Carbon in 
27 plates of GaAs made from a single rod were 
analyzed by the dry method in different machine 
times of two cyclotrons. The mean results were 
24.6 ppb wt with cr=2.S ppb. For the last five 
plates taken from adjacent parts of the rod, the re- 
suits were 23.25±0.84 ppb. We separate "F by pre- 
cipitation as KB'*F4 [2]. This method has shown 
satisfactory precision. 



2.3 Calibration of IR Spectrophotometry 

Our calibration curve for IR spectrophotometry 
of C in Si is shown in figure 1 [1]. About 70 care- 
fully prepared samples were submitted to round- 
robin IR measurement in 24 organizations. About 
1/3 of them were then analyzed by CPAA; some of 
them were also analyzed by SIMS, The differences 
between the results of CPAA and SIMS are proba- 
bly caused by the presence of background in SIMS. 
Satisfactory calibration curves were obtained also 
for O and N in Si [3,4]. Efforts are now being made 
to obtain a reliable calibration curve for C in GaAs. 



3. Accuracy for ERDA and NRA for H 

Figure 2 shows ERDA spectra of a H- and D- 
containing amorphous silicon film on a Si wafer 
[5], The detection efficiency is seen to depend 
highly on the detection angle. No suitable sub- 
stances are found as comparators for H and D, and 
their quantities are usually obtained in comparison 
with an internal standard element by the use of the 
theoretical relationships and experimental parame- 
ters. The H content of amorphous silicon films ob- 
tained by ERDA and vacuum fusions are given in 
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